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1 Introduction 


The most significant hazard associated with spacing of commercial aircraft 
when landing is the wing tip trailing vortex from the leading aircraft. The 
vorticity associated with the circulation that generates lift on a wing is con- 
tinually shed and tends to roll up into a tight vortex aligned with the flight 
path. There are two of these vortices for a given aircraft, one associated 
with each wing. The air velocity fluctuations associated with these vortices 
are strong enough to pose a threat to following airplanes. Thus, the FA A 
regulations for spacing of aircraft at landing are intended to account for 
the decay of these vortices. This decay, however, is not well understood. 
As described by Stough et al. [1], the current FAA regulations depend on 
whether the prevailing conditions are VMC (visual meteorological condi- 
tions) or IMC (instrument meteorological conditions). Under VMC, it is up 
to the pilot of the following aircraft to maintain a safe distance, whereas 
FAA minimum spacing requirements must be adhered to under IFC. Thus, 
the only way in which the meteorological conditions are taken into account 
is in regard to the pilot’s visual capabilities, not the dynamics of the trailing 
vortices. It is believed that with more knowledge of how trailing vortices 
decay under varying atmospheric conditions, air traffic capacity and safety 
may be increased. 

Past research has primarily focused on how the trailing vortices move un- 
der stable or neutral atmospheric conditions. Greene’s [2] analytical model 
compared well with the few laboratory and observational measurements 
available, but dealt only with wake vortices in the free atmosphere, not 
in the atmospheric boundary layer. Zheng and Ash [3] have had success 
with their numerical model predicting the motion of the vortices near the 
ground. In their case, however, only the turbulence associated with the vor- 
tices, themselves was considered. To date, there has not been a successful 
study describing the interaction of realistic atmospheric convective turbu- 
lence and the trailing vortex system. The goal of our research is to fill this 
gap- 

We will use the TASS(Terminal Area Simulation System) model devel- 
oped by Proctor [4] to study this problem. TASS is a non-hydrostatic Large 
Eddy Simulation model which includes parameterizations for hydrometeors. 
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2 Objectives and goals of the research 

The proposed research involves four tasks. The first of these is to simulate 
accurately the turbulent processes in the atmospheric boundary layer. TASS 
was originally developed to study meso-7 scale phenomena, such as tomadic 
storms, microbursts and windshear effects in terminal areas. Simulation of 
wake vortex evolution, however, will rely on appropriate representation of 
the physical processes in the surface layer and mixed layer. This involves two 
parts. First, a specified heat flux boundary condition must be implemented 
at the surface. Using this boundary condition, simulation results will be 
compared to experimental data and to other model results for validation. At 
this point, any necessary changes to the model will be implemented. Next, 
a surface energy budget parameterization will be added to the model. This 
will enable calculation of the surface fluxes by accounting for the radiative 
heat transfer to and from the ground and heat loss to the soil rather than 
simple specification of the fluxes. 

The second task involves running TASS with prescribed wake vortices in 
the initial condition. The vortex models will be supplied by NASA Langley 
Research Center. Sensitivity tests will be performed on different meteorolog- 
ical environments in the atmospheric boundary layer, which include stable, 
neutral, and unstable stratifications, calm and severe wind conditions, and 
dry and wet conditions. Vortex strength may be varied as well. Relevant 
non-dimensional parameters will include the following: Richardson number 
or Froude number, Bowen ratio, and height to length scale ratios. The model 
output will be analyzed and visualized to better understand the transport, 
decay, and growth rates of the wake vortices. 

The third task involves running simulations using observed data. MIT 
Lincoln Labs is currently planning field experiments at the Memphis airport 
to measure both meteorological conditions and wake vortex characteristics. 
Once this data becomes available, it can be used to validate the model for 
vortex behavior under different atmospheric conditions. The fourth task will 
be to simulate the wake in a more realistic environment covering a wider 
area. This will involve grid nesting, since high resolution will be required in 
the wake region but a larger total domain will be used. 

During the first allocation year, most of the first task will be accom- 
plished. 
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3 Work accomplished during the period 1/94- 
6/94 

3.1 Surface energy balance scheme 

The parameterizations to be used in the surface energy balance have been 
determined. A description of these parameterizations follows. 

The surface energy balance may be described as: 

A Q, = Q*-Q 9 -Qh-Qe (1) 

where A Q 8 represents the rate of change in energy storage in the surface 
layer of the soil per unit area, Q* the net radiative flux to the ground, Q g 
the heat loss to the deep soil, Q h the sensible heat flux to the atmosphere, 
and Q e the latent heat flux to the atmosphere. 

The radiative flux may be broken down in the following way: 

Q* = K + Ll-L T (2) 


where L j is the longwave incoming radiation and L | is the longwave 
outgoing radiation. K is the net shortwave radiation and may be written as 

K = TkS{ 1 - A) sin $ (3) 


where S is the solar radiation at the top of the atmosphere, A is the albedo, 
and Tk is the atmospheric transmissivity. ® is the solar elevation angle and 
may be expressed as 


(4) 


In this formulation, <(> is the latitude, A c the longitude, 6 S the solar declina- 
tion, and tuTC the standard time of day at 0° longitude. 

Burridge fc Gadd [5] parameterize the transmissivity as 


T k = (0.6 + 0.2sin$)(l - 0.4 (tc h )(1 - 0.7<r Cw )(l - 0.4 <t C l ) (5) 


where <^c H , a C l represent high, middle, and low cloud fraction, 

respectively. 

We may express the outgoing longwave radiation as 

L T= €<tT* (6) 


3 



where e is the emissivity of the ground, a the Stefan-Boltzmann constant, 
and T g the ground temperature. 

Staley & Jurica [6] parameterized the incoming longwave radiation as 

L i= [<t c + (1 - <t c )(0.67)(1670 9o ) 008 ]<tT o 4 (7) 

where <r c is the total cloud cover, q a is the humidity at the surface refer- 
ence level of the atmosphere (usually about ten meters), and T a is the air 
temperature at that level. 

By Monin-Oboukov scaling, the sensible heat flux to the atmosphere is 


Qh = pc v u m 6+ 


( 8 ) 


where 


ku a 

Um ~ 


( 9 ) 


and 


„ k(0a ~ * g ) 


( 10 ) 


Here k represents Von Karmann’s constant, Zo the roughness length, 
L the Monin-Obukhov length (which depends upon stability), and a-nd 
tl>H are stability dependent functions related to the dimensionless wind shear 
and potential temperature gradient. u a and 6 a are the velocity and potential 
temperature, respectively, at the reference level in the surface layer, and 9 g 
is the potential temperature at the ground. 

There are two remaining terms on the right hand side of Equation 1, 
Q g and Qe* To determine these, we will use a slab model of the soil which 
was developed by Bhumralkar [7]. Deardorff [8] showed that this was a very 
accurate method. In addition, it is much more efficient than numerically 
solving the heat conduction equation with many layers within the soil. As 
shown in Figure 1, the soil is modeled as two layers: a slab in which the 
temperature changes throughout the day, and a substrate in which the tem- 
perature remains constant throughout the diurnal period. The amount of 
heat energy in the ground slab lost to the substrate is 


Q ' = 2*p,d 1 '. {T '_ Tm y (11) 

Here, p 9 is the soil density, c M is the soil heat capacity, d\ is the depth of the 
ground slab, T g is the temperature of the ground slab, T m is the temperature 
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Figure 1: Schematic drawing of slab soil model. 


of the substrate, and ti is the diurnal period (twenty-four hours). The rate 
of change of the energy contained within the ground slab (the left hand side 
of Equation 1) is then 

A Q, = PsCsdi^- (12) 

In a manner similar to Equation 8, we may define the latent heat flux 
to the atmosphere as 

Qe = —pL e u m q+ (13) 

where L c is the latent heat of evaporation, 


_ %a ~ ft) 

9 * {Mff)- 


(14) 


and q a and q g are the specific humidity at the surface and ground level, 
respectively. Using Deardorff’s [8] parameterization, 


9a — <]g — a/ [9a ” 9«at(^») 1 


where q 8a t(T a ) is the saturation specific humidity at the ground temperature 
and 


a = min(l, — ). 

Wk 


(16) 
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Here, w, is the volume fraction of soil moisture in the ground slab and ttffc is 
the soil moisture fraction above which the ground acts as if it were saturated. 
We may then write 




(17) 


where E g is the evaporation, P the precipitation, and w m the moisture 
fraction in the substrate. C\ and C 2 are empirical constants. 

After we are convinced that TASS gives the proper results with a given 
surface heat and moisture flux, this surface energy parameterization will be 
tested with data from Ripley & Redmann [9], which includes measurements 
of u a , <5*, Q g , Qe, and Qh- We will prescribe the shortwave radiative flux, 
the surface wind, the stability, moisture, and the deep substrate temperature 
to match the experiment and follow the development of T ty Qe , and Qh- 
We will then add this parameterization to the TASS model. 


3.2 Large eddy simulations of the atmospheric boundary 
layer 

Since May 18, 1994, Dr. David Schowalter and Mr. David DeCroix have 
been in residence at NASA Langley Research Center learning how to use 
the TASS model under the direction of Dr. Fred Proctor. They will remain 
there until mid-August. What follows has been accomplished there in the 
first month. 

In order to validate model results for simulation of the atmospheric 
boundary layer, a surface heat flux boundary condition has been added to 
the model. The surface heat flux in this case is a function of the time of day. 
The model has been tested with input data from the Wangara field experi- 
ment [10]. We are using data from Day 33 to test mean quantities against 
experimental values. The surface heat flux for this case was not measured 
and was obtained from Deardorff’s [11] model. Only instantaneous values 
of velocity and temperature are available from this observational data set 
for altitudes above 16m. Thus, we are comparing turbulent quantity results 
with Deardorff [11] and other modelers who tested large eddy simulation 
output against the same experimental data. 

Figure 2 shows a comparison of potential temperature profiles as a func- 
tion of the time of day. The simulation was initialized with the profile shown 
at 9:00, which matched the experimental data exactly. A 40x40x40 grid was 
used for the internal domain. Periodic horizontal boundary conditions were 
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Figure 2: Variation of potential temperature throughout the day. Compar- 
ison is with data from the Wangara Experiment, Day 33. 


used with a rigid lid and sponge condition at the top boundary. The domain 
covered a region bKm x 5 Km X 2Km, the latter dimension representing the 
domain height. This corresponds to 125m resolution horizontally and 50m 
vertically. This resolution corresponds to that used by Deardorff, which 
was desirable for comparison purposes. The simulated profiles at 12:00 and 
15:00 compare well with the data. It is important to note, however, that 
our results show slightly higher temperatures at the top of the mixed layer 
and a larger mixed layer height than those measured at 15:00. DeardorfPs 
LES model results show similar discrepancies. This is most probably due to 
large scale advection and subsidence taking place during the observations. 
These synoptic effects are not accounted for in the current model, though it 
may be in the future if deemed necessary. 

Figure 3 reveals contours of vertical velocity in a vertical plane at 14:00 
(2:00 p.m.) for the Wangara simulation. Thermal plumes can be seen quite 
clearly. It should be mentioned at this point that the size of visible turbulent 
eddies is limited on the small scale end by the resolution of the grid (125 m 
in the horizontal). The thermals shown here have a typical horizontal size 
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Figure 3: Vertical plane showing contours of vertical velocity w at t = 14 : 00 
for the Wangara simulation. Contours are from -1.4 m/s to 2.2 m/s by 0.2. 
Negative contours are shown with dashed curves. 
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Figure 4: Horizontal plane showing contours of vertical velocity in at t — 
14 : 00 and at a height of z — 199m. Contours are from -2.0 m/s to 3.25 
m/s by 0.25. Negative contours are shown with dashed curves. 
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Figure 5: Same as Figure 4, but showing horizontal velocity vectors at 
z = 224m. 


10 




Figure 6: Horizontal plane showing contours of vertical velocity w at t — 
14 : 00 and at a height of z = 1149m. Contours are from -1.5 m/s to 4.25 
m/s by 0.25. Negative contours are shown with dashed curves. 
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Figure 7: Same as Figure 6, but showing horizontal velocity vectors at 
x = 1174m. 
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Figure 8: Sensible heat flux profiles as a function of local time of day. (a) 
Current results. (6) Results from Deardorff [11]. 
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of roughly three grid cells. Thus, it is possible that the characteristic size 
of eddies would be smaller if a smaller grid size were used. These types of 
sensitivity tests will be done in the future. These thermals can be observed 
from above, as well, as in Figure 4, which shows a horizontal plane of vertical 
velocity contours. Figure 5 shows the horizontal velocity vectors at nearly 
the same height. Careful observation will show that at this altitude, the 
thermals are associated with horizontal convergence. Figures 6 and 7 show 
similar contours and vectors, but higher up in the boundary layer (1149m 
and 1174m, respectively) where the thermals are less intense. At this higher 
altitude, there is divergence associated with the thermal plumes. Funda- 
mentally, then, the simulation gives a realistic representation of boundary 
layer processes. 

In figure 8, we compare the sensible heat flux profiles (not measured 
in the experiment) with DeardorfFs results. Qualitatively, they are very 
similar, though the heat flux minimum in these curves has a slightly lower 
value in our case. Because we are comparing model to model, however, 
it is very difficult to determine which results are correct. It is important, 
therefore, to compare turbulent quantities with observational data. 

4 Work in progress and objectives for the period 
7/94-10/94 

Because we believe that the wake vortex development will be sensitive to at- 
mospheric turbulence, it is important to test turbulent quantities in addition 
to mean values. The Wangara data does not include these measurements, 
so it is necessary to turn to the Minnesota field experiment [12]. Most im- 
portantly, we will compare variances of velocity and temperature as well as 
turbulent fluxes of heat, moisture, and momentum between the model and 
the experiment. Not only is this a useful test of the TASS model specifically, 
but of large eddy simulation in general. To our knowledge, no one has done 
an exhaustive comparison between atmospheric turbulence data and large 
eddy simulation results in the past. 

After these tests are performed, we will add the surface energy balance 
to the model so that the surface fluxes will be determined by the amount 
of solar radiation reaching the ground and by the heat transfer through 
the soil. This work will be under way by the end of the allocation period 
(October 31, 1994). The first task is expected to be completed entirely in 
March, 1995. 
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